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single leptonic Jarlskog invariant as predicted in the absence of new physics effects. 
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1 Introduction 



In recent years the experimental study of neutrino oscillations has much contributed to 
our knowledge of particle physics by establishing non vanishing neutrino masses and by 
measuring or constraining the corresponding mixing angles. In the near future, within 
the domain of neutrino oscillations, the most important tasks for experiments, beyond 
refining the measured values that are already known, are the measurement of the mixing 
angle ^13, which at present is only limited by an upper bound, the determination of the 
sign of Amg]^ = mg — ml, which would solve the existing ambiguity between normal 
(positive sign) and inverse (negative sign) hierarchy in the neutrino mass spectrum, and 
finally the observation of CP violation. Once a positive signal for CP violation in neutrino 
oscillations has been established, the next question will be, like in the analogous case for 
quarks, whether the observed phenomena are in agreement with the simplest picture where 
all CP violating quantities in oscillation phenomena are described in terms of only one 
phase 6 appearing in the PMNS mixing matrix, or equivalently in terms of the Jarlskog 
invariant J, which in the standard notation for mixing angles is given by: 

J = — cos sin 26'i3 sin 26'23 sin 26'i2 sin 5. (1) 
8 

As well known, in the mass matrix of 3 Majorana neutrinos two additional CP violating 
phases also appear. As they do not contribute to CP violation in neutrino oscillations 
we will not consider them here, but they can have observable effects in other phenomena 
like neutrino-less double beta decay and leptogenesis. A lot of theoretical effort has been 
devoted to develop and discuss the strategies for measuring 5 in long baseline experiments 
with Super-Beams, (3- beams and Neutrino Factories. The most suitable observables for 
this purpose are the asymmetries associated with the "golden" channel z/g ^ or Uf. v^. 
However other channels are also important to disentangle a set of discrete ambiguities 
("degeneracies") in the solutions for ^13 and 5 that arise in the procedure: the "intrinsic" 
degeneracy due to the functional dependence of the asymmetry on 6*13 and 5, the "sign" 
degeneracy that occurs if the sign of Amgj^ is unknown and the "octant" degeneracy that 
depends on the sign of ^^23 — 7r/4. Each ambiguity is 2- fold, so that a total of 8 solutions 
are generated if all options are open. The intrinsic degeneracy would remain even if the 
signs of Am|]^ and of 6'23 — 7r/4 will be known at the time of the measurement. The next 
important channel, the "silver" channel v^. ^ or z/g ^ Vj. can be very useful to remove 
the ambiguity, as well as the i/^ 01 z/^ ^ Vr- In this respect also the possibility of 

varying the L/ Ey ratio of the baseline lenght L and the beam energy E^, has been discussed 
in the literature. 

In the present article we address the related problem of testing the PMNS mechanism 
for CP violation (indicated here as the "standard" framework). We want to study the most 
suitable tests and discuss the required sensitivity. In order to get a quantitative estimate 
of the possible effects of new physics we will consider the new physics contributions to 
the relevant observables in a model where the PMNS unitarity is relaxed. This model, 
denoted as the Minimal Unitarity Violation model (MUV), has been introduced in refs. 
[T] and further studied in [2], where the present bounds on the new parameters of the 
model have also been derived. Although this model is certainly not general it is however 



sufficiently structured to provide an indicative estimate of the possible deviations from the 
standard description of leptonic CP violation that one may expect0. Different models were 
studied in the literature with new physics interactions at the source and at the detection 
of neutrinos in oscillation experiments [5] and also during the propagation in matter [6]. 
All these models refer to new interactions, with both diagonal and non diagonal lepton 
flavour properties, parameterized in terms of effective 4-fermion vertices with couplings 
left free but constrained by existing data. In the MUV model there are less parameters, so 
that there are more rigid constraints from the data. The resulting MUV framework is one 
that reproduces the broad features of the more general models but where the parameters 
at the source, at detection and in matter are related. As the role of MUV for us is to 
show that, after imposing the existing experimental bounds, there is still a lot of space 
for new physics to affect the values of the CP violating observables, the fact that there 
is even more freedom is not a problem. But the presence of new physics can also have 
additional implications. For example we will assume in the following that, by the time 
of the CP violation measurements, some quantities will be known within a reasonable 
precision, the most important one being the value of 6'i3. Then, in the presence of new 
physics, the results of these measurements could be distorted by the new effects so that 
they do reproduce the correct values of the relevant quantities. For this purpose, even when 
the model independence of the results is guaranteed by MUV, it could be invalidated in a 
more general model. We will discuss this aspect in the following. 

The conclusion of our investigation is that the same set of measurements that have 
been discussed for the determination of 6 and for the removal of degeneracies offer a viable 
and effective testing ground for the mechanism of CP violation. As already mentioned, we 
assume that the value of ^13 will be known with sufficient precision in a model independent 
way by the time of the measurement of CP violation and, as a consequence, we work with 
^13 fixed at a set of representative values (sin ^13 = 0.02, 0.05, 0.10, 0.15). Within the 
MUV model the golden channel probabilities are the least affected by new physics, given 
the powerful constraints arising from the present and forthcoming limits on the fi —>■ ej 
decay, and can be used for an essentially model independent determination of sin and 
also of 6. If the measured asymmetry is compatible with a discrete set of values of 6 (if not 
compatible one directly obtains clear evidence for some exotic mechanism) then the silver 
channel and other possible measurements offer the possibility to check the corresponding 
predictions of the standard picture. For discussing this comparison we eliminate S and 
express one CP violating quantity in terms of another. For example, a CP violating 
asymmetry for the silver channel can be given in terms of the measured value of the 
golden channel asymmetry, or the value of one asymmetry can be predicted in terms of 
the same asymmetry at a different L/Ej^ value. We show that the present constraints 
on the MUV parameters still allow a considerable space for results at sizable variance 
with respect to the standard expectation, so that meaningful tests can be obtained with 
realistic measurements. For this to be realised the experimental uncertainty associated 
with the CP observables in the standard model must be smaller than the deviations induced 

^Previous studies of new sources of CP violation in the context of non- unitary neutrino mixing matrix 
can be found in 3 and 



by the new physics parameters. To check that this is indeed the case, we estimate the 
errors on the various asymmetries at one particular experimental setup (a neutrino factory 
with L = 1500 Km) showing that, even under reasonable assumptions on efficiencies, 
backgrounds and systematic errors, one can reasonably expect to incisively probe the 
presence of new physics. We also show that the degeneracies do not prevent the possibility 
of testing the standard CP mechanism, at least for not too small values of ^13. 

This paper is organised as follows. In Sect. 2 we briefly review the standard formalism. 
In Sect. 3 we summarise the MUV model and extend the formulae for probabilities and 
asymmetries to include the new physics effects. In Sect. 4 we discuss the set of possible 
future facilities that we will consider in our analysis. Sect. 5 contains our main results on 
the way of testing the standard picture of CP violation. Finally in Sect. 6 we summarise 
our conclusion. 



2 Transition probabilities and asymmetries in the PMNS 
framework 

In this section we briefly recollect the formalism used to describe standard neutrino os- 
cillations in matter and in vacuum. The lepton mixing is described by the 3x3 unitary 
matrix Upmns |Z] which, in its standard parameterization, has the following form [8], [9j: 



PMNS 



C12C13 S12C13 Si3e~'^ \ 

-C23S12 - Si3S23Ci2e" C23C12 - SisS23Sl2e^^ S23C13 
V S23S12 - Si3C23Ci2e^'^ -S23C12 - Si3C23Si2e'^'^ C23C13 J 



(2) 



where Sij and Cij are short-hand notations for sin6'ij and cos%, respectively whereas 6 is 
the (Dirac) CP violating phas^. Neutrino oscillation probabilities are obtained by solving 
the Schrodinger equation for flavour eigenstates with the effective hamiltonian given by: 

H = ^UpMNsdmg{0, Aml^, Aml^)Ul,MNS + diag(\^, 0, 0) (3) 

where Amfj = mf — m|, and diag(y, 0, 0) is the potential felt by neutrinos in their propa- 
gation in matter. V is given by: 

V = V2GfN^ 

with Gp being the Fermi constant and A^^e the ambient electron number density [lOj. 
Expressions of the transition probabilities P^i'a ^p) = Pap {.P{.^a ^p) = Pi 



apj 



m 



matter can be found in a number of papers (see, e.g., [TT]-[l3]). All our numerical results 
in Sect. 5 have been obtained by numerically solving the coupled Schroedinger equations 
involving the appropriate Hamiltonians (for the SM and MUV cases) . On the other hand, 
for the sake of simplicity, the analytical formulae, given in the present section as well as 



^Possible Majorana phases are not explicitly considered here because they do not manifest themselves 
in oscillation experiments. 



in the next one, have been computed using some useful approximations, motivated by the 
underlying physics. 

At present neutrino oscillation results for the mixing angles [H], within l-cr, are con- 
sistent with the so-called tri — bimaximal mixing (TBM) matrix |15j 



^ PAINS 



/ V3 ^ \ 

1_ J_ J_ 

\ J- L ±. I 

\ V6 V3 V2 / 



(4) 



Following [16], we parametrize the deviation from tri-bimaximal mixing by introducing 
three small parameters r,s,a in such a way that: 

r 1 1 

5l3 = ^, si2 = -^(l + s), S23 = ^(l + a). (5) 

The results of global fits on neutrino data ([I7],[18]) imply that these parameters can be 
at most of O (10%). One can then expand all probabilities and asymmetries in r,s,a 
around the TBM matrix. We can further simplify the analytical expressions by expanding 
in the variable ^ = Am2i/| Am|^| ~ 0.03 and working in the small A21 limit, where 
A21 = Aml-^^L/AEi, <^ 1. In the following we quote the relevant transition probabilities 
and the related CP asymmetries up to second order in those small parameters. Note that 
only r appears in the asymmetries up to second order while a and s appear among the 
dominant terms only in P^^ (eq. [181 ESI [27]) and in P^^ (eq. [26l) . 
First, we consider the golden channel Pe^ |12j : 



, sin^ (A - 1) A31 4 ^ ,^ ^ X . / . ,x ^ sinfAAai) 

4 2 sin^(AA3i) 
9^ 

where A = VL/2A^i. The CP conjugate channel Pgp is obtained from Pg^t by changing 
S —6 and A —A. The resulting expression of the asymmetry Aefj, in matter is quite 
cumbersome; we then limit ourselves to explicitly quote only the vacuum case which shares 
some of the relevant features of the corresponding quantity in matter: 

^ _ 12r A21 sin^ sin^ A31 

~ 4Aii + 9r2 sin2A3i + 6rA2i cos(5 sin2A3i 

Notice that the "T-conjugate" asymmetry A^^ can be easily calculated by assuming 
that the density profiles are symmetric under interchange of the position of the neutrino 
source and detector, so that Pa/3 {S, A) = P^a (— ^, A)^ which in turn implies A^e {^1 ^) = 
Aefj, (—5, A). Thus, the results for P^e and A^e can be directly obtained from eqs.([6l[Zl)- 

In the approximations adopted in this section, the silver channel Per [19] has an ex- 
pression similar to eq.Q, the only difference being a sign change in the second term where 
the CP violating phase appears. This implies: 

^ _ 12r A21 sin^ sin^ A31 

~ -4Aii - 9r2 sin^ A31 + 6r A21 cos5 sin2 A31 ^' 



Finally, we quote 



Af,r = -rAaisin^. (9) 

This simple expression is a consequence of the zeroth-order contribution of P^r (the ex- 
pression is given is App jA.3|) . which is absent in Pe^,T- 

Notice that all the asymmetries in vacuum are suppressed by the small quantities 
A21 and 9i3; however, since the denominators in eqs.([7}lH]) are also suppressed, a partial 
cancellation is at work and, in particular, one generically expects Ag^, A^t ^ A^^-. 

Matter effects play an important role in modifying the behaviour of the asymmetries 
as function of the CP phase S (see, e.g., [20] and [21]). In fact, while in the vacuum 
case Aap = when sin 5 = 0, this is no longer true in the more general case because 
the passage through matter introduces fake CP-violating effects in neutrino propagation. 
Then, to extract genuine CP violating effects, one often defines the subtracted asymmetries 
as A^^^{6) = Aa/3{6) — Aa/3{6 = 0). We prefer to deal with more directly measurable 
quantities, so that we use unsubtracted asymmetries which, for non negligible matter 
effects, are non vanishing at 5 = 0. 

To conclude this section, we want to mention that also the CP-conserving channel P^^ 
can be in principle used to extract information on the CP phas^ because of the term 
r A21 cos(5sin2A3i, which shows the same suppression as A^r- In particular, in [22] it has 
been shown that the variation of P^^ due to a change in the value of 6 can be as large as 
30% when L > 2000 Km and < 2 GeV. 



3 Summary of the MUV model 

As a representative framework for physics beyond the standard model affecting neutrino 
oscillations we consider the Minimal Unitarity Violation model (MUV) introduced in [T] 
and further studied in [2]. This framework is based on the possibility that, although the 
complete theory is unitary, the effective leptonic mixing matrix Upmns describing neutrino 
oscillations at low energy is not unitary, as is the case in a large class of theories addressing 
the question of neutrino masses. The deviation from unitarity is parametrized in terms of 
a Hermitian matrix r] in such a way that (using the same notation as in [2]): 

= (1 + r^) UpMNS 

where A^ is the non-unitary leptonic mixing matrix relating the interaction and mass 
eigenstates, Ua = Nai Vi- Being a hermitian matrix, ri contains 9 new parameters, for 
example 6 moduli and three phases. In particular, the phases are new sources of CP 
violation in the leptonic sector. The moduli of the ri matrix elements are bounded at the 
level of (9(1%) (or smaller) from universality tests, rare lepton decays, the invisible width 
of the Z boson and neutrino oscillation data [1]: 

We\ ~ I^mmI ~ l^rrl < 2.5 X 10"^ 



^See app |A.4l for details. 



|r/ej < 3.6 X 10^^ 

\r]er\ < 8 X 10"^ 

\t]^^\ < 4.9x10-1 (10) 

On the other hand the phases, as for the standard CP phase in Upmnsi are completely 
unbounded. As mentioned in the Introduction, in models where new physics contributions 
are introduced at all stages in the oscillation process, starting from the source of neutrinos 
and then in the interaction with matter during the propagation and finally in the detection 
process, there is space for a wider range of effects. In fact the new physics allowed by MUV 
is more limited than in general. But our main purpose in showing the MUV predictions 
is to prove that, in each specific case, there can be sufficiently large new physics effects to 
be detected. On the other hand, while we have specified the set of present constraints, it 
must be kept in mind that more restrictive limits could be obtained by the time when the 
CP violation experiments will be performed. 

3.1 Transition probabilities and asymmetries in MUV 

In the MUV case, since our aim is to propose consistency checks of the leptonic CP sector 
of the standard model based on the currently planned experimental options, we present 
analytic transition probabilities and asymmetries with the same approximations used in 
SectEl Furthermore, we also expand up to second order in the small r^a/s's even if, from 
the limits given in eq.f fTOi) . the retained rj terms could be smaller than the neglected terms 
in the r, s and a expansion (defined in eq.(l5])). The asymmetries including the MUV 
corrections are derived using the procedure described in Appendix A of [2]. Based on 
the discussion of the previous section, in the following we only quote the vacuum results, 
which are enough to understand the role of the MUV parameters compared to the standard 
model results. 

In the MUV framework, the probability Pe^ mainly depends on the parameters rje^ and 
?7eT- and their related phases. The detailed expression is given in App JA.l] here we only 
retain terms which are linear in the small Tj'^pS and do not depend on the product rj^p A21, 
obtaining: 



Pe^ = P^^ + 7]e^,r sin A31 [3 sin(5 - A31 - ^e/.) + sin(5 + A31 - 4^)] 

+2rieTr cos(5 — der) sin^ A31 (11) 

where by P^ we denote the SM results. The asymmetry A^^ obtained from eq. ffTTl) reads: 

(121 



A 



12 r A21 sin 5 sin^ A31 + 18 r r/e^ sin 2A31 sin(5 — 5, 



4 A21 + Qr^ sin^ A31 + 18r sin^ A31 [rjer cos((5 — Ser) — Vefi cos{6 — Se^)] + 6r A21 cos 6 sin 2 A31 



Similarly to the standard model case, the modified can be obtained from Ae^ changing 
the sign of all the phases. We see that a new physics term proportional to r/e^t appears 
in the numerator with a different dependence on L/E^\ however, in the MUV model, due 



to the strong bounds on the parameter rjefj,, we can safely neglect this term, containing 
a new CP violation effect. Note also that the new terms in the denominator show a 
dependence on L/Ei, which is the same as the standard terms so that their contribution 
(and in particular that due to r/er) cannot compete against r unless is very small. As a 
result the asymmetry does not deviate much from the standard result. It follows that, 
in the MUV model, this channel is particularly suitable for an experimental determination 
of the UpMNS parameters, as we now briefly discuss. But it must be kept in mind that in 
a more general model there can certainly be new physics effects in leading to larger 
deviations from the standard picture of CP violation than in MUV. 

In the following analysis we assume that the value of will be known with sufficient 
precision by the time that a meaningful measurement of the CP violating phase 5 can 
be performed. As we are contemplating the possibility of important new physics effects, 
an obvious question is whether 6'i3 can be determined in a reasonably model independent 
way. In the MUV framework, Pg^ and the corresponding asymmetry A^.^ are optimal in 
this respect as the new physics effects are already bound to be reasonably small. The 
planned experimentation at T2K |23] indeed plans to measure 6'i3 through this channel. 
Alternatively, reactor experiments like Double Chooz [2l], would measure ^13 from the 
diagonal transition probability Pgg which, in the limit of vanishing rj^.^, is given by: 

Pee = l-2sin2A31 

In this channel the new physics effects can indeed be large. Thus a measurement of 6'i3 
from reactors only would not guarantee a model independent determination even in the 
MUV case. But the combination of reactors with T2K would answer the question on 
model independence and actually could put bounds on the maximal value allowed for rjer 
in the MUV model. In a more general context than MUV different measurements are even 
more necessary in order to put limits on the model dependence of the measured value of 
6'i3. In the following we will assume that the value of ^13 is known, but it is well possible 
that the issue of a truly model independent determination of ^13 will only be solved by 
the same generation of experiments that will tackle the CP violation problem. 

Important corrections to the SM results may appear in Per which, for rje^ = 0, is given 

by: 

Per = P!^ + Ver {r siu A31 [3 sin(5 - 6er " A31) + sin(5 - 6er + A31)] - 

^ A31 e [sin(2A3i - 6er) - 3 sin 6er]] + i vl (5 + 3 cos 2A31) (14) 

The new source of CP violation is only linearly suppressed in the small rjer parameter; 
in particular, the first term in curly brackets can, in principle, compete with the standard 
term cos((5 — A31) A31 sin A31 because it is not suppressed by the small quantity A21. 
Note also that the quadratic term corresponds to the zero-distance effect discussed in [1] . 
From Per, the asymmetry Aer reads: 



2 I 2 



+ 2rrier cos(5 — 6e 



Q ^1 



(13) 



(15) 



A^ 



-8 r A21 sin 6 sin^ A31 + 4 rjer [3r sin 2A31 sin(5 — 6er) + A21 sin 5er (3 + cos 2A31)] 



8/3 A|i + 6 r2 sin^ A31 - 4r]er A21 sin 2A31 + 3r cos(5 - 6er) sin^ A31 - 4r A21 cos 6 sin 2A31 



where, for simplicity, we have disregarded terms proportional to rjl^ and rjeT^2i- It is clear 
that, whenever A21 ~ r^er, the asymmetry can be quite different with respect to the SM 
result. Then A^t can be a good CP- violating observable to perform consistency checks of 
the leptonic CP sector of the SM. 

A quite different case is that related to A^^; as illustrated in eq.®, in the SM A^^- is 
somewhat limited to small values, around 10"'^ for the maximum allowed and 5. Thus 
it would be much more easy to see any deviation from standard results if new physics 
would affect transitions [2j. This is indeed the case since the main correction to 

the SM expression of P^r is linear in the new parameter 77^,-: 

P^r = P^r' - 2 V>.r sm 2 A,, siu V (16) 

and, consequently: 

A^^r = - 4 cot A31 sin 5^r (17) 

where the SM term is suppressed by at least one power of A21. The new physics term 
easily overwhelms the standard one and A^,- can reach values as large as O (10~^). 

Let us finally comment on P^^ in the MUV scheme. The correction to the SM results 
is quadratic in the small parameters of our perturbative expansions and only depends on 
1]^,^ (but not in 77^^): 

^MM = ^/f/f + 4 Vt^r cos S^r sin^ A31 (2 O + 7]^^ COS V) (18) 

4 A list of possible future experiments 

In this section we summarize the experimental facilities we consider in our analysis. In re- 
cent years there has been a vast amount of theoretical work on the performances of several 
setups aimed at studying some of the still unknown neutrino parameters and properties, 
like the value of ^^13, the existence of leptonic CP- violation, the type of mass hierarchy 
and the precise value of 623- In particular, to assess the problem of CP violation, it has 
been pointed out that the zeroth-order requirement is to have very intense neutrino fluxes 
and good detector performances: this automatically aims at second-generation neutrino 
experiments, like Super-Beams, /3-beams and Neutrino Factories. Since the measurement 
of the standard CP phase 6 is directly connected with the solution of the problem of the 
degeneracies [25]- [28], one needs to combine different transition channels (golden, silver 
and also z/^ —>■ Vr channels), or to exploit several L/E^j regimes (which can be accom- 
plished either using two detectors located at different baselines or by varying the neutrino 
energy). Also, a combination of these facilities among themselves and/or in association 
with atmospheric neutrino data can be quite useful for constraining 5. For a review of 
these possibilities, see [29]. For our purposes, it is sufficient to select some of the most 
promising future experiments for the measurement of 5 and probe their ability to perform 
the consistency checks we aim to. The extension of the analysis to other facilities and 
combinations is straightforward. 
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Figure 1: The standard model asymmetries Aefj,, A^t for different values of sin 613 = S13 = 
0.02, 0.05, 0.10, 0.15 as a function of the phase 6. We have assumed that the sign of 
A31 is positive (normal hierarchy), that 623 is maximal and sin^ 6*12 = 1/3. We have also 
fixed the mass differences to Am2i = 8 x 10~^ eV^ and Am^i = 2.4 x 10"'^ eV^. Neutrino 
energy and baseline distance L are fixed to E = 30 GeV and L = 1500 Km, respectively. 

• Super-beams 

These are conventional neutrino beams from vr decay (except for their large intensity). 
Since the flux is predominantly i/^, the and transition channels 

are available (the Vf^^ transition can only be probed using high energy neutrino 
fluxes having sufficient energy to produce r leptons). The second generation of 
upgraded Super-Beams can run with both neutrino helicities, thus having sensitivity 
to CP violation; we refer to two different options: the upgraded version of T2K [23j, 
T2HK, with L = 295 Km and {E^) ~ 0.75 GeV and the CERN Super-Beam project 
SPL |30] with {E^) ~ 0.3 and L = 130 (Frejus) or L = 732 (Gran Sasso); 

• /3-Beams 

These are intense z/e and z/g fluxes produced from boosted radioactive ion decays [31]: 
then it is possible to study z/g — z/g, i^e — ^ i^fi (and probably z/e Vr if very high- 
energy options will be availabe) and their CP-conjugate channels. Depending on the 
choice of decaying ions and boost 7 factors, several possibilities could be exploited. 
Here we refer to the standard low-energy /3-beam, with i^Ej) ~ 0.4 GeV and L = 130 
Km (LE/3B) 02] and high-energy p-heam (HE/3B) with {E^) ~ 1 GeV and L ~ 732 
Km |32]-[33l. 




0.001 




-0.004 ' ' ' ' ' ' ^ ' ' ' ' ' ^ 

2 4 6 2 4 6 

5 

Figure 2: The standard model asymmetry A^r for different values of sin 613 = S13 = 
0.02, 0.05, 0.10, 0.15 as a function of the phase 6. All other parameters are fixed as in 
the previous figure. Neutrino energy and baseline distance L are fixed to E = 35 GeV and 
L = 1500 Km, respectively. 

• V Factories 

In a neutrino factory, neutrinos are produced via the decay of muons, which are 
accumulated in a storage ring and accelerated to the desired energy. The result- 
ing neutrino flux is analitically calculable from three-body decays kinematics and 
strongly depends on the parent muon energy. Since two different types of neutrinos 
are produced for fixed muon charge, there are in principle six different channels that 
can be simultaneously probed, assuming that the appropriate detector technologies 
are available. Then, in principle, the three asymmetries Ag^, and can be 
investigated, as well as the disappearance. If we concentrate on the sensitivity 
to CP violation, the performance on a neutrino factory especially depends on the 
values of ^13. For smaller values, the main uncertainties in the 5 extraction come 
from the degeneracy problem and a good performance is reached for i?^ > 30 GeV 
and L G [3000 — 5000] Km [34] (especially if in combination with another experiment 
running at the "magic" baseline L ~ 7500 Km [35]) whereas for larger values matter 
effects are the main source of uncertainty and a smaller baseline is preferabl^. Here 
we consider a neutrino factory with = 50 GeV, for which {E^^^p^) = 30 GeV and 
= 35 GeV. 

''Notice also that a low-energy neutrino factory [36] (i?^ = 4.12 GeV) seems to have a good performance 
to CP violation for intermediate 6*13 ~ 3°. 
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Table 1: Summary of the relevant channels and parameters of the future neutrino facilities 
considered in this work. In the second column we list all the possible available channels at 
a given facility, assuming that the adequate detector technology will be available. 

In TablH we summarize the setups considered in our analysis. 
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Figure 3: Scatter plots, for the absolute values of the asymmetries Ae^ and A^f, of the 
domains spanned when the MUV parameters are moved in their allowed ranges. The angle 
^13 is fixed to sin^^i3 = 0.1. The values of the other parameters are as in the previous 
figures. Energies and baselines match the values quoted in Tab)^ As explained in the 
text, in the MUV model the deviations from new physics are particularly small for these 
asymmetries. 
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5 Testing the CP violation mechanism 

We start by recalling the dependence of the A^^, A^t and A^^ asymmetries on the PMNS 
phase S in one typical long baseline experiment with L = 1500 Km and E^, = 30 GeV 
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Figure 4: Scatter plots, for the absolute values of the asymmetry A^r of the domains 
spanned when the MUV parameters are moved in their allowed ranges. The angle 613 is 
fixed to sm6'i3 = 0.1. The values of the other parameters are as in the previous figures. 
For any of the two facilities, we consider = 35 GeV. 
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Figure 5: Scatter plots, for the absolute values of the asymmetry A^r of the domains 
spanned when the MUV parameters are moved in their allowed ranges. The angle 613 is 
fixed to sin6'i3 = 0.1. The values of the other parameters are as in the previous figures. 
For any of the two facilities, we consider Ey = 30 GeV. 



(for Ae^ and Aer) or E^ = 35 GeV (for A^r)- In Figs {T] and [2] we plot the asymmetries vs 
6 for fixed values of 6*13, sin 6*13 = 0.02, 0.05, 0.10, 0.15, assuming that the sign of A31 



is positive (normal hierarchy), that ^23 is maximal and sin^ 6'i2 = 1/3. We also fixed the 
mass differences to Amg^ = 8 x 10~^ eV^ and Amg^ = 2.4 x 10~^ eV^. The curves refer 
to the conventional picture (all MUV parameters set to vanish). We see that, at fixed 
L/ Ey, for a given value of each asymmetry there is in general a two-fold ambiguity in the 
corresponding value of 5. 
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Figure 6: The asymmetry A^t plotted as a function of A^.^ for sin^is = 
0.02, 0.05, 0.10, 0.15. We have assumed that the sign of A31 is positive (normal hi- 
erarchy), that 023 is maximal and sin^ 6^12 = 1/3. We have also fixed the mass differences 
to Am2i = 8 X 10~^ eV^ and Amg^ = 2.4 x 10~^ eV^. In each panel we use L = 1500 
Km and Ey = 30 GeV. Notice the different x,y-axes scale in each plot. In each panel 
the marked closed line describes the standard model relation between the two asymmetries, 
while the dots span the domain obtained by varying the MUV parameters. 

In Figsl3],|4]and[5]we show the domain spanned when the MUV parameters are moved in 
their allowed ranges, for each of the three asymmetries (in absolute value) and sin ^13 = 0.1. 
The asymmetry can be measured at both /3-beams and Neutrino Factories whereas 
the Super-Beams are better suited to measure A^e; in particular, we illustrate the effect 
of considering the MUV parameters at NF@1500, T2HK, HE/3B and SPL@732, which 
have quite different values of L/{ES) (from L/{ES) ~ 50 Km/GeV for NF@1500 up to 
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Figure 7: The asymmetry A^r plotted as a function of A^n for sin ^13 = 
0.02, 0.05, 0.10, 0.15. The values of the other parameters are as in the previous fig- 
ures. We compute A^^ for E^, = 30 GeV and A^j. for = 35 GeV, both at L = 1500 
Km. Notice the different x,y-axes scale in each plot. In each panel the marked closed line 
describes the standard model relation between the two asymmetries, while the dots span the 
domain obtained by varying the MUV parameters. 



L/{E^) ~ 2440 Km/GeV for SPL@732; notice also that T2HK, SPL@130 and LE/?B share 
almost the same Lj {E^) ~ C(300 — 400)). We can easily see that the prediction of Aefj, at 
each 6 value is not much affected by varying the MUV parameters in their experimental 
allowed range at any of the facilities that we considered. 

This is not the case for the asymmetries Afj_r and Aer'- while in the standard picture they 
are bounded in absolute value, these bounds can be grossly violated in the MUV model. 
In particular, the measurement of A^^j. can in principle be done at Neutrino Factories with 
different L/{E^)] as we can see in FigJU the strongest effects of new physics are obtained 
at the facility with the smallest L/{Ey), thus showing that the Neutrino Factories with 
small baselines and large neutrino energies are better suited to detect new physics effects 
in the —>■ Vr channel [2]. 

Finally, in Figj5]we present the results for the asymmetry AeT- A part from the matter 
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Figure 8: The asymmetry Ag^- plotted as a function of A^.^ for sinOi^ = 0.10, for a Neutrino 
Factory with L = 1500 Km. The values of the other parameters are as in the previous 
figures. The solid line refers to the standard model prediction, with the uncertainties for 
both the asymmetries evaluated at some representative values of the phase S. Dots represent 
the spread of MUV predictions for the same quantities. 

effects wliicfi modify tfie sliape of tfie asymmetry, there is not mucli difference between 
tlie two scenarios, as a consequence of tlie fact tliat tlie values of L/{Ey) differ at mucli 
by a factor of 3. 

Clearly if tfie experimental values of one or more asymmetries are not compatible with 
the standard bounds, the PMNS model is directly disproved. Instead we consider here the 
case where the measured values are compatible with the bounds. From the measured value 
of 5 the prediction for all CP violating asymmetries is obtained in the standard model. It 
is then interesting to eliminate 6 and plot each asymmetry as a function of the measured 
value of Ae^, which is the least affected by the MUV new physics. The resulting plot is in 
the form of a closed line. If the standard model is valid the measured values for each pair 
of asymmetries should fall on the line. These plots are shown in the four panels in Figj6] 
for Af,T vs Af.^ and in FigJ7]for A^^. vs at L = 1500 Km and = 30 GeV (for A^,^ 
and Aer) or E^, = 35 GeV (for A^^), always for sin6'i3 = 0.02, 0.05, 0.10, 0.15, a positive 
sign of A31, 623 maximal and sin^6'i2 = 1/3. These plots show that even with moderate 
precision a meaningful test of the PMNS mechanism of CP violation can be achieved. In 
fact the spread which is still allowed by the existing bounds on the MUV parameters is 
quite large for A^t and A^^. In more general frameworks for new physics the deviations 
can even be larger. In particular the interval of allowed values for A^fj, could be wider. 
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Figure 9: The asymmetry A^j. plotted as a function of A^.^ for sin 613 = 0. 10, for a Neutrino 
Factory with L = 1500 Km. The values of the other parameters are as in the previous 
figures. The solid line refers to the standard model prediction, with the uncertainties for 
both the asymmetries evaluated at some representative values of the phase 6. Dots represent 
the spread of MUV predictions for the same quantities. 



5.1 A discussion of the expected experimental accuracy 

A central question now is related to the experimental capability to get a good enough 
measurement of the asymmetries at the facilities that we consider: if the CP violating 
quantities will not be measured with sufficient precision, then the test can be invalidated 
because one cannot distinguish the deviation from the standard model results due to MUV 
effects. To illustrate this point, we consider a neutrino factory with L = 1500 Km and 
En = 50 GeV and build the relevant integrated asymmetries from the number of expected 
events Nf^ and iVg in a given detector according to: 



A 



where the event rates for the Va ^p and the CP conjugate Ua 
computed from: 

Nfs= [ dE, P^p{E,) ap{E,) ^{E,) ep{E, 



dE,P^-p{E,) (TpiE.j — 



■{El,) ei3{Ey 



(19) 

VjB transitions are 
(20) 
(21) 



in which o"^(/3) is the cross section for producing the lepton /3(/9), Sp^p) the detector efficiency 
to reveal that lepton and (t)a{a) the initial neutrino flux at the source. The asymmetries 
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Figure 10: Effect of the octant and sign degeneracies on the A^^ and A^r asymmetries, for 
sin ^13 = 0.02, 0.05, 0.10, 0.15. In each panel we used L = 4000 Km and = 30 GeV. 
The solid line refers to 6*23 = 42° and normal hierarchy, the dashed one to 6*23 = 42° and 
inverted hierarchy, the long-dashed to 623 = 48° and normal hierarchy and crosses refer to 
023 = 48° and inverted hierarchy. 



we are interested in are Ae^, A^t and A^^; then f3 can be a muon or a tau (and their 
antiparticles) . We perform a simulation of a neutrino factory in the spirit of [12], taking 
1 ■ 10^^ /i"*" X 4 years, 1 ■ 10^^ fi^ x 4 years [37] and putting all the detectors we need 
at the same baseline L = 1500 Km. It is obvious that different technologies are needed 
to detect muons and taus; for this reason, we imagine to look for golden muons using an 
improved version of the MIND detector described in [38] and, to evaluate the uncertainty 
associated with A^^, we use the efficiencies and backgrounds from |39j. The detector mass 
is fixed to 50 Kton. The search of r events takes place looking for its n decay mode; in the 
following we use an Opera-like detector (the so-called Emulsion Cloud Chamber (ECC) 
detector With this detector, the silver channel Ue has been carefully studied 

in [19], and we refer to these papers for our estimate of efficiencies and backgrounds. 
However, such a detailed study is still missing for the z/^ — > Ur transition in the context 
of the neutrino factory; here we follow [2] and assume that it should also be possible to 
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Figure 11: New physics effects in the MUV model in presence of the sign degeneracy for 
Aefi and Aer asymmetries, computed for sin 9i3 = 0.1 at L = 4000 Km and E^, = 30 GeV. 
Solid lines refer to the standard model scenario whereas dots are the MUV predictions. 

look for the other r decay modes, thus gaining a factor 5 in statistics (at the prize of 
increasing also the backgrounds) □. The mass of the detector is 10 Kton. An important 
point to stress is that at the level of precision required to see new physics effects, systematic 
uncertainties play a relevant role; in the following we assume a conservative overall 10% 
systematic error for the silver detector and an overall 2% for the golden one. For the 
sake of illustration, in the following we integrate the rates in eqs.(l20l)-(!2Tl) just over the 
whole range; this means that the asymmetries do not profit of any spectral feature of 
the signal (certainly available with such good-performance detectors), which would help 
in extracting more accurate information on CP- violation. In this respect, our results can 
be improved in a more refined analysis which is beyond the scope of this section. The 
results are shown in FiglHl where we plot the asymmetry as a function of and 
their corresponding standard model uncertainties, at some representative value of the CP 
phase 6, as well as the spread predicted in the MUV environment. The value of 613 is fixed 
to sin 6*13 = 0.10. It can be clearly seen that the asymmetry uncertainties are not as big as 
to spoil the possibility to see substantial deviations from the standard model predictions, 
since many of the values obtained in the MUV framework are well beyond the asymmetry 
uncertainties For comparison, this scenario is equivalent to the third panel of FiglBl 
however, the numerical values of the asymmetries are different, as a result of integrating 
the transition probabilities over different neutrino and antineutrino cross sections which, 
at the relevant neutrino energies, are in the ratio crp/cr^ ~ 0(2-3). 

similar option has been also investigated for the silver channel first in 34J and further discussed in 

^Notice that the uncertainty on A^r is larger than that on j4e^, as a consequence of the smaller event 
rates and larger background for the silver channel, although this is not clearly visible in the figure because 
of the different horizontal and vertical scales used to make the plot more understandable. 
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Figure 12: Effect of the octant and sign degeneracies on the A^^ and A^j. asymmetries, for 
sin^is = 0.02, 0.05, 0.10, 0.15. A^e is computed for L = 295 Km and = 0.75 GeV and 
Aer is computed for L = 1500 Km and E^ = 30 GeV. The solid line refers to 623 = 42° 
and normal hierarchy, the dashed one to 623 = 42° and inverted hierarchy, the long-dashed 
to 023 = 48° and normal hierarchy and crosses refer to 623 = 48° and inverted hierarchy. 



The same can also be done with the other asymmetry, A^.^. This is illustrated in 
FigJHl in which we plot A^r as a function of Ae^ and their corrresponding uncertainties. 
It is clear that this case looks less favourable, as while the size of the experimental errors 
is comparable the MUV scatter plot is more concentrated. Still, also in this plot, there 
are MUV predictions for A^^ not falling inside the error band (whereas, as expected, 
the already small deviations for are completely overshadowed by the uncertainty). 
However, compared to FigJTl A^r is now substancially different from zero; in particular, 
the almost constant value is a direct consequence of the ap-to-a^ ratio TZ; in fact, assuming 
for simplicity that Nf^ TIP and iVg ~ P, we can estimate A^r as follows: 

^ (^-^)p+_(p-p) (7^-l) 
Tzp + p ^ (n + i) 



where in the last step we have assumed that P — P ~ in the standard model. Thus, 
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Figure 13: New physics effects in the MUV model in presence of the sign degeneracy for 
A^e, computed for L = 295 Km and Ey = 0.75 GeV, and A^t asymmetries, computed for 
L = 4000 Km and Ey = 30 GeV, at fixed sin6'i3 = 0.1. Solid lines refer to the standard 
model scenario whereas dots are the MUV predictions. 

IZ ~ 0(2-3) justifies the almost constant value of in FiglHl We stress that our 
error estimates can be improved in several ways; in particular, we have checked that the 
uncertainties on and A^^ are rather sensitive to the assumed systematic errors and this 
is especially important if one want to amplify the MUV effects in the Uj. transition. 

As a general conclusion, this exercise shows that the indications that can be drawn by 
looking at the asymmetries at the probability level, computed at the peak energy of the 
neutrino spectra, are not invalidated when reliable error estimates are taken into account. 
It is also true that, a part from the ECC detector, the proposed facilities have not been 
built yet and MonteCarlo simulations for the detector response are still in progress. For 
these reasons, we prefer to investigate the other relevant issues of the paper using the 
non-integrated asymmetries, as done at the beginning of this section. 

5.2 The impact of degeneracies 

One may wonder what would be the impact of the sign and octant degeneracies on the 
significance of the tests if the corresponding ambiguities were not already fixed at the 
time of the CP violation experiment. This is studied in Figs. [TOlfTTl In Fig{TO]we see the 
possible impact of the sign and the octant degeneracies for the asymmetries computed in 
the standard picture for a very long baseline experiment with L = 4000 Km (so that the 
matter effect is particularly large and the sign degeneracy is amplified) and Ey = 30 GeV. 
In this plot for the octant degeneracy we have assumed that 623 = 42°. The results show 
that the ambiguity induced by an octant change is not very important because the 623 
angle is in any case close to maximal. The sign ambiguity instead makes a lot of difference 



at such a long baseline distance. This implies that an experiment of this type can indeed 
easily measure the A31 sign if this is not already known. In fact, as can be seen from FigITT] 
(where, for simplicity, the curves are plotted only for ^23 = 42°) the new physics effects 
would not be large enough to fill up the gap arising from the sign change, at least for not 
too small values of 6^13 . Conversely, the sign degeneracy does not spoil the possibility of a 
meaningful test for the presence of new physics. 
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Figure 14: The asymmetry Ae^ computed at two different baselines L = 4000 Km (hori- 
zontal axis) and L = 1500 Km (vertical axis) at the same neutrino energy E = 30 GeV, 
for sin = 0.02, 0.05, 0.10, 0.15. As for the other parameters A31 is positive, 623 is 
maximal and sin^^i2 = 1/3. 

For a shorter baseline the sign degeneracy can induce a larger confusion. This is 
illustrated in Figs JT2lfT3l that correspond to Figs JTOlfTTl but with A^e computed for L = 295 
Km and E^, = 0.75 GeV and A^r computed for L = 1500 Km and E^ = 30 GeV. In this 
case, first of all, a better precision would be needed to disentangle the sign ambiguity. 
But we see that, barring the unfortunate case that the new physics effects are just right 
to shift the data point from one closed curve to the other one, a meaningful test of the 
standard picture is still possible. 

Another complementary possibility of searching for new physics effects that could also. 
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Figure 15: The asymmetry A^r computed at two different baselines L = 4000 Km (hori- 
zontal axis) and L = 1500 Km (vertical axis) at the same neutrino energy E = 30 GeV, 
for sin = 0.02, 0.05, 0.10, 0.15. As for the other parameters A31 is positive, 623 is 
maximal and sin^^i2 = 1/3. 

to some extent, bypass the sign degeneracy problem is the comparison of a given asym- 
metry at different basehne distances. This comparison would be possible in principle 
even if only the golden channel could be measured. The plots in Fig{T3] show A^^ mea- 
sured at L = 4000 Km and = 30 GeV versus the same asymmetry at L = 1500 Km 
and El, = 30 GeV for sin6'i3 = 0.02, 0.05, 0.10, 0.15, A31 positive, 6*23 maximal and 
sin^ 6'i2 = 1/3. We see in a clear way that, for this asymmetry in the MUV case, the 
possible new physics effects are not sufficiently pronounced to allow their detection in an 
experiment with realistic precision. This confirms that while A^^ is the best suited asym- 
metry for a model independent determination of the phase 6, for a test of the standard 
model one needs to also measure one or more additional asymmetries. The situation is 
completely different for the Aer asymmetry measured at the two indicated baselines, as 
shown in Fig{T5] obtained for the same values of all parameters. Here the effects of new 
physics are large enough to be detectable with reasonable experimental precision. 



6 Conclusion 

In this paper we have considered the proposed facihties for detecting CP violation in 

neutrino oscillation experiments. All of these experiments aim at measuring the value 
of S. We have studied the next question of testing whether the observed CP violation is 
consistent with the standard model prediction that all CP violation observables in neutrino 
oscillations must be proportional to the leptonic Jarlskog invariant. As a quantitative 
model of new physics we have adopted the MUV framework with the present bounds on its 
parameters. Although the MUV model is a rather restricted model of new physics (larger 
effects can be expected in a more general model) , we have found that indeed the deviations 
from the standard picture of CP violation could be detected with realistic experimental 
accuracies (namely those accuracies which are needed to make a significant measurement 
of S and to disentangle the related sign degeneracy). Actually our results show that the 
deviations allowed by the present constraints on the MUV parameters are large enough 
that a meaningful test can be done with a relatively modest precision. We have explicitly 
checked this statement by evaluating the uncertainties associated with the asymmetries 
Af,i^, AgT- and A^^ as if they were measured at a neutrino factory with detectors at L = 1500 
Km, showing that new physics effects can produce deviations from the standard model 
asymmetries well outside of the error bars (in this respect A^t is much more promising 
than Af^T- as the asymmetry to be measured in correspondence with Ag^j). But it is also 
true that the data constraining the new physics parameters will improve in the near future, 
so that by the time that CP violation can be observed the allowed new physics deviations 
will presumably be more contained. For testing the PMNS mechanism of CP violation 
it is necessary to measure not only the golden channel A^.^ (which is the most protected 
from the new physics effects in the MUV framework) but at least one more channel, for 
example the silver channel A^t (or A^,-). One can either measure both channels at the 
same value of L/E^ or the silver channel at two (or more) values ol L/E^,. In the first 
relatively small value of L/E^, is to be preferred, because the ratio of the effects 
from new vs standard physics in the asymmetries increase with decreasing L/E^,. Thus 
L/Ei, must be sufficiently large to allow a good counting rate for asymmetries, so that 
CP violation can be observed, and sufficiently small to enhance the visibility of the new 
physics effects. From our simulations it appears that a good compromise corresponds to 
a neutrino factory with L = 1500 Km and E^j^ ~ 50 GeV. 
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A Transition probabilities in the MUV framework 

In this appendix we quote the vacuum probabilities obtained in the MUV framework (and 
in the standard model if not previously quoted in the main text) using the approximations 



described in Sectsj2] and [31 



The transition probability Pg^ reads: 

Pe^ = + Ve^ {r sin A31 [3 sin((5 - A31 - 5^^) + sin(5 + A31 - 5^^)] 

+-A21 [sin(2A3i - 4^) - 3 sin 4^]} + 2r]er [r cos(5 - 4^) sin^ A31 (22) 
2 

+-A21 cos(A3i - 6er) sin A31 } + 

r^e^r^er sin A31 [-3 sin(A3i - 5er + Se^) + sin(A3i - Se^, + Ser)] + 
IvU^ + 3 COS 2A31) + vl sin' A31. 

Let us analyze the new CP violating terms, namely those linearly suppressed by the 
new parameters rje^ and rjer (we can safely assume that the second-order term rje^rjer 
is smaller than the others). In each curly brackets there is a competition between the 
term suppressed by 6^13 and that suppressed by A21; since A21 ~ (9(10~^) L/ E^, for small 
enough L/Ey it can be neglected, unless r is of the same order of magnitude (in fact, this 
is almost the regime in which we work because L/ E^ is at most (9(10^) Km/GeV). For the 
coefficient of rjer this means that we are left with a CP conserving term and it does not 
contribute to CP violation in the e/i channel. On the other hand, for r/g/^, the relevant term 
is CP violating and only linearly suppressed with ^13. We recover here the observation 
made in [2] that at short enough baselines, the new CP violating terms in P^^ mainly 
depend on the phase associated to the corresponding rjap- This also means that, due to 
the strong bounds on ?7e^ in MUV, is a good CP violating observable to measure the 
standard phase 5. 



A.2 P, 



er 



The MUV effects in the Per probability are obtained from Pg^ of eq. (lA.ip changing r^g^ ^ 

Tier and (5e^ ^ 5 

Per = Pif + Ver {r siu A31 [3 sin((5 - A31 - 6er) + Sm{6 + A31 - 5er)] 
2 To 

+ -A21 [sin(2A3i - 6er) - 3 sin 4^]} + 2r]e^ |r cos(5 - ^g^) sin A31 (23) 
2 

-I--A21 cos(A3i - 5e^) sin A31 } + 

r]ef,r]eT sin A31 [-3 sin(A3i + 6er - Se^,) + sin(A3i + 5g^ - 6er)] + 
1 
2' 



-rilib + 3 cos 2A31) + ril^ sin^ A31. 



We can apply to Pgr the same considerations as for Pg^, concluding that the relevant 
new CP violating effects are all contained in the terms proportional to rjer- 



A. 3 P, 



The probability reads: 



P^, = P^^ - 2ri^r sin 5^r { sin 2A31 - ^ A21 cos 2A31 



r cos(5 - 6er) sin^ A31 + - A21 cos(A3i + 6er) sin A31 



+ 



(24) 



r]f,r sin 6^r sin 2A3i(?7^^ + r]rr) 
where the standard model expression is given by: 



P^^ = sin^Asi + -A2isinA3i(rsin(5sinA3i + SCOSA31) +sin2A3i(4a2 -r^) (25) 

-^A2i sin2A3i - ^A^i(l -3cos2A3i). 

Contrary to Pe^ and Per, the larger new physics effect in P^,- is only linearly suppressed 
by r]^r- This allows to approximate the probability as in eq. fll6p and to get the simple 
expression for the asymmetry A^r as given in eq. (fT71) . 



A.4 P, 



The expression of the transition probability P„„ reads: 



MM 



-Pu/f + 4 VfiT cos 6^r sin^ A31 (2 a + r]^r cos 6f,r ] 



(26) 



where 



P 



SM 
MM 



1 - (1 - 4a^) sin^ A31 - - A^^ (1 + 3cos2A3i) + ^ ^21 sin2A3i (1 - s - r cos 5). 

(27) 
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